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We analyze the Hertz-Moriya-Millis theory of an antiferromagnetic quantum critical point, in the 
marginal case of two dimensions {d = 2, z — 2). Up to next-to- leading order in the number of 
components (A^) of the field, we find that logarithmic corrections do not lead to an enhancement of 
the Landau damping. This is in agreement with a renormalization-group analysis, for arbitrary A''. 
Hence, the logarithmic effects are unable to account for the behavior reportedly observed in inelastic 
neutron scattering experiments on CeCue-ajAu^. We also examine the extended dynamical mean- 
field treatment (local approximation) of this theory, and find that only subdominant corrections to 
the Landau damping are obtained within this approximation, in contrast to recent claims. 

I. INTRODUCTION 

Materials in the vicinity of a quantum critical point (QCP) continue to be the subject of intensive investigationgi. 
Particular attention has been paid to the heavy fermion system CeCue-^Aua; which undergoes a phase transition 
from a paramagnetic, heavy-fermion metal to an antiferromagnetic metal as a function of chemical composition and 
pressure^. Another remarkable case is the stoichiometric compound YbRh2Si2 which is apparently located very close 
(on the magnetic side) of an antiferromagnetic QCB^. Hydrostatic pressure stabilizes the magnetic phase in that case. 
As the critical point is approached, these systems exhibit an enhancement of the specific heat coefficient C /T, magnetic 
susceptibility, resistivity Ap and thermoelectric power, which is not easily described within the conventional theory 
of a three dimensional quantum critical antiferromagnet. 

Neutron scattering experiments^i^i^ on CeCug-ajAuj; have revealed that the spin fluctuation spectrum of this 
compound has a two-dimensional character in a wide range of temperatures. Indeed, maxima of the neutron scattering 
intensity are observed along rod-like structures in reciprocal space along which spin fluctuations are independent of 
one component of the wavevector. This motivated Rosch et al. « to propose that the observed anomalies are the results 
of the coupling of three-dimensional electrons to quasi two-dimensional critical spin fluctuations. In such a situation, 
conventional Hertz-Moriya-Millis theorj&Siifliii of quantum critical behavior leads to C/T cx ln(l/T), Ap cx T, 
consistent with experimental findings. A similar observation applies to the thermoelectric power^^. The microscopic 
reason for the two-dimensional nature of spin fluctuations is not clear however for CeCug-xAua;, while magnetic 
frustration might provide a natural explanation^ in the case of YbRh2Si2 (for which inelastic neutron scattering data 
are not yet available). 

The approach of Rosch et al^ requires flne tuning on a microscopic level (since a three- dimensional dispersion 
of the electrons will generically lead to three- dimensional critical fluctuations at the QCP) but it is internally 
consistent on a theoretical level. Indeed, the vertex corrections to the electron- spin fluctuation interaction is flnitei^ 
in the case of a coupling to a three-dimensional fermion spectrum. (Note that the situation is different for two- 
^ dimensional fermions, which leads to a singular vertesti^). Since the vertex is non-singular, the physics of long- 
k> I wavelength magnetic fluctuations is described by a (j)^ fleld theory with ohmic damping. The quartic term is marginally 
5—1 ' irrelevant for d — 2, which is the upper critical dimension of this theory. Hence, the flow of the quartic term will lead 
^ to logarithmic deviations from mean-field critical behavior. These effects have not been investigated previously for 
the frequency-dependent response function in this context. Inelastic neutron scattering2*Sil on CeCue-^Aua; at the 
QCP {xc — 0.1) have revealed an anomalous enhancement of the Landau damping. Whether logarithmic effects at the 
upper critical dimension could account for such an enhancement in the range of experimentally accessible frequencies 
is an outstanding open question which we address in this paper. We use the number of components (N) of the field 
as a control parameter, and perform a calculation of the dynamical susceptibility to order We find that, to this 

order, the sign of the corrections to Landau damping cannot explain the experimental data. 

An alternative theoretical viewpoint on the physics at the QCP is that the spin-fluctuation self-energy at the 
QCP is purely local (i.e momentum independent) and has an anomalous power-law dependence on frequency and 
temperature. This was first proposed2i& as a phenomenological fit to the inelastic neutron scattering data in the form: 
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IrnXglq = uj°'f{Lu/T) with a ~ 0.7 — 0.8. Si et ali^ have further developed this point of view in the framework of 
the extended dynamical mean-field (EDMFT) theory of the Kondo lattice^^'iS'ii. It has been argued in this context 
that the separation of electronic degrees of freedom and long-wavelength magnetic fluctuations is not legitimatai^^, 
because the (Kondo) coherence scale below which electron degrees of freedom can be eliminated might vanish at the 
QCP. Recently however, Grempel and Si showediS that starting above the Kondo temperature the EDMFT of the 
Kondo lattice can be mapped, using bosonization tricks, onto the EDMFT of the (j)'^ theory for magnetic modes (see 
also the earlier work by Sengupta and one of us^). They further suggested that, in two dimensions, this theory 
would lead to an anomalous Landau damping, with the power-law form mentioned above {a < 1). We study this 
effective theory in the EDMFT approximation, using both a strict 1/N expansion and a numerical solution using 
the large-N expansion as an "impurity solver". We find that, contrary to the statements of Ref. |19|, no anomalous 
power-law is generated. In fact, the EDMFT approximation to the 0"* field-theory in the marginal case does not 
correctly reproduce the logarithmic corrections to the damping rate found in our direct calculation for d = 2, and 
leads only to subdominant 0(u;^) corrections at low-energy. This extends to the marginal case the previous work of 
Motome and two of us^*' on the the EDMFT approximation applied to critical behavior, in which it was indeed found 
that the approximation is better in higher dimensions, and becomes unreliable below the upper critical dimensions. 

Hence, the general conclusion of our paper is that logarithmic effects at the upper critical dimensions in the con- 
ventional theory of an antiferromagnetic QCP cannot explain the behavior observed experimentally for CeCug-a^Au^:, 
at least at dominant order in the 1/N expansion. 

This paper is organized as follows: in section ^ we perform a general renormalization-group analysis of the 0** 
theory with damping. In section ITTll we take a large- iV perspective, and present an explicit solution up to next to 
leading order in the 1/N expansion. We also describe the local (EDMFT) approximation to this theory, up to order 
In section Hvl we discuss in more details the EDMFT approach, and solve numerically the EDMFT equations 
using a large-N impurity solver. We conclude with a discussion of the implications of our work for the understanding 
of non-Fermi liquid behavior in heavy-fermion materials close to an antiferromagnetic QCP. 



II. RENORMALIZATION GROUP ANALYSIS 



A. Model and general renormalization-group framework 



In this paper, we consider the following field theorySiiSiii for an A^-component field (pa (a = 1, • ' ' ; -^)- 

'^ = ^E^o'('?'*^)E'^'('?'*^) + IT / d'xdr[Y,M^,r)y (1) 

q.uj a a 

in which the free propagator reads, for imaginary frequencies: DQ^{q,iuj) = r + \lu\ + q^. This theory describes the 
vicinity of an antiferromagnetic quantum critical point, corresponding to the (bare) value z = 2 of the dynamical 
critical exponent. In this context, the dissipative term is induced by the coupling to electronic degrees of freedom 
(Landau damping). The vector is the difference from the ordering wavevector, and r is a parameter which measures 
the distance to the critical point. We are particularly interested in this paper in the two-dimensional case, but we 
shall also briefly consider lower dimensions 1 < d < 2 (Sec. lII(!)l . In this section, we shall use the following normalized 
coupling constant: 

50 ^ S,uo with : S, = ^r(rf/2)(4.)./2 (2) 

When considering the bare theory, an (ultra-violet) cutoff A on momenta is introduced. We consider the bare 
irreducible 2-point function F = which is the inverse of the correlation function (dynamical susceptibility 

x{q,iuj) = {(j>{q,iu!)(/){q, —iuj))). The renormalization group (RG) specifies^^ how F changes upon a rescaling of 
the cutoff A —^ AA. The corresponding RG equation is obtained by imposing that the renormalized 2-point function 
= ZT is independent of the original cutoff (with Z the field- or "wave function" renormalization). Focusing first, 
for simplicity, on the field-theory in the critical (massless) case, F satisfies the following equation (all Green's 
functions in this and the next subsection are at T = 0): 



T[q,iLo;gQ,A]^{XAfZ^'T 



1 ■ A 1 

AA'(W'"'^'' 



(3) 



In this expression, Z\ and the running coupling constant gx are defined from the usual RG functions: 



A^5A = /3(gA) , 5A=i=5o (4) 
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•^37 ^^^^ ^(^a) . ^A=l = 1 



(5) 



Let us emphasize a key aspect of Eq. Q , namely that the damping term does not require an independent renormal- 
ization so that the frequency dependence in the r.h.s of Q involves only the wave-function renormalization Z\. This 
can be proven to hold to all orders, using a field theoretic RG^^. For a use of this method in a different context see 
ReflU 

We have calculated the RG function /3 at one-loop order and the function 77 at two- loop order. The diagrams 
contributing to the 2-point function F up to this order are depicted in Fig. ^ We obtain (with e = — 2): 

m = -eg + &2.g' + 0{g^) 
v{g) = C2.g' + 0(5^) 

where the coefficients 62 and C2 read, in the N-componcnt theory: 



b. = 



N + 8 



C2 = 



6 ' 144 
Details of calculation of the coefficient C2 are given in Appendix ^ 



(6) 
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L 




FIG. 1: Graphs contributing to the 2-point function at two-loop order. 



B. The marginal case d = 2: logarithmic corrections to Landau damping 

We focus here on the marginal case c? = 2 at the quantum-critical point, and integrate the above RG equations in 
order to obtain the corrections to the Landau damping at T = 0. Integrating the RG equations 14I5|I using © yields, 
to this order: 

90 

9\ 



1 - 62.90 In A 



Z), = exp ( \ (7) 

\ 1 — 62.90 In A / 

As expected in the marginal case, gx flows logarithmically to zero as A — > 0, while Zx tends to some (non-universal) 
constant Z*. To analyze the frequency dependence of F we use the general RG equation setting q = and 
choosing A = A* such that lo Z x" / {X* A)'^ ~ 1- This leads to: 

r{q = 0,iu;go,A) = \uj\^{gx^) (8) 

in which we use the notation '^{g) = T{q — 0,iuj — i; g, A = 1). We then expand the r.h.s of (O in powers of gx*, with: 
^A' ^ —1/(62 In A*) ^ — 2/(62 ln(Z*|aj|/A^)]. This expansion actually starts at second-order in gx, since the 2-point 
function F is subtracted to insure that one sits at the critical point (in other words, the tadpoles Ii do not contribute 
to the frequency dependence). Noting that: 12(9 = 0,iuj = 0;gx,A) — 12(9 = 0,iuj,gx,A) = cg1\uj\, coefficient c 
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computed in Appendix ^ we finally obtain a correction to the frequency dependence of the 2-point function at the 
T = QCP of the form: 

vt n • ^ \ [67r^ln2-llC(3)] (TV + 2) 1 ^ 

Hence, only subdominant corrections to the Landau damping are generated in the marginal case d = 2. Furthermore, 
the positive coefficient in the above expression can be interpreted as an effective exponent a — 2/z > 1. This is in 
agreement with the 1/N expansion, which we consider later in this paper. 



C. e = 2 — d expansion of critical exponents and corrections to scaling 

We consider here the case d < 2, in which the quartic term is relevant, and the effective coupling gx tends to a 
non-trivial fixed point. The case d = 1 is relevant, for example, when consideringii an Ising chain in a transverse field 
(TV = 1) or a chain of quantum rotors (N > 1) coupled to a dissipative environment. 

We make use of the general RG equation (jSJ by choosing A = A* such that g/A*A = 1. For A ^ 0, we now 
have: gx — > g* and Zx ~ A** (with the critical exponent rj = ri{g*). Hence, we obtain (at T = 0) for the dynamical 
susceptibility x — a-t low-frequency and small momentum: 

xiq,uj) = q-'+''H^/cq^-'') (10) 

In this expression, is a universal scaling function associated with the fixed point. From this expression, the dynamical 
critical exponent is identified as: 

z = 2-'n (11) 

An expression which holds to all orders, and is the result of the existence of a unique RG function, as explained above. 
This scaling form can actually be generalized to finite temperature in the quantum critical regime in the form: 



with $ a universal scaling function, and ci a non-universal constant. In particular, at zero-momentum, Tx{q = 0,'^) 
is an entirely universal scaling function of tu/T. At small lo/T, we expect an analytic dependence on lu/T, with 
T{0,uj)/T = Ci — iC2Uj/T + . . .. In contrast, at large uj/T we have r(0,ti;)/T — —iC^uj/T, and the subdominant 
terms are not analytic, but are related to those in H19|l below. Here Ci_3 are all non-trivial universal constants. The 
value of C3 is determined in the e expansion by the computations in Appendix \K\ (C3 = 1-1- c(e/52)^ -I- . . .), and this 
universality is clearly related to the universal logarithmic correction in The values of C1.2 require a computation 
of the damping at T > 0, and accurate determination of these is likely to require a self-consistent treatment of the 
loop corrections, as discussed in Ref'23. In particular, C2 ^ C3 and hence 

lin, lim ^J^^£(^ ^ lin, lin, ^J^£i^_ (13) 

So the Landau damping co-efficient is sensitive to the order of limits of a; — > and T — > 0; we expect that similar 
considerations will also apply to the logarithmic corrections in d = 2 noted in J^. 

From the RG expressions in Section III Al the e-expansion of the critical exponent rj = 2 — z reads, to order e^: 

(A^ + 2)(12-7r2) , , 
We have also obtained the correlation length exponent at order e: 

Finally, we comment briefly on the next-to-leading corrections at T to the Landau damping at the QCP. These 
are obtained from corrections to scaling, i.e. depend on the manner in which gx flows to the fixed point: 

9x^9* + i9o~g*)X'' + --- (16) 
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with: 

n = p'{g*)=e + 0{e^) (17) 
The RG equation then yields the following form for the 2-point function, including corrections to scaling: 

r(g,a;) = g2-, [^^{^q--) ^ q^^^^{u:q-^) ^ . . .] (18) 

Taking the limit g ^ 0, this yields the following form of the Landau damping term: 

T{q^Q,iLo)^Co\Lof^ +C2\uj\^^ +■■■ (19) 

In the present case, the exact identity z = 2 — r/ implies that the dominant term is simply cx \lo\, while a correction 
of the form |a;| "^2-,, jg found. The logarithmic correction found above can be seen as the limiting behavior of this 
correction in the marginal case d — 2 (e = 0) . 



III. LARGE-N ANALYSIS OF THE (P^ THEORY OF AN ANTIFERROMAGNETIC QUANTUM 

CRITICAL POINT IN TWO DIMENSIONS (2 = 2) 

In this section, we treat the field-theory within the large-N expansion. For this purpose, we scale the coupling 
constant in by 1/A^ and set: 

Uo = ^ (20) 

We denote by Eq and Ei the contributions to the self-energy of order and 1/A^, respectively {{(j)(t))~^ = Dq^ — T,). 
The (skeleton) diagrams contributing to Eq and Ei are depicted in Fig. |21 The corresponding expressions read: 

J2j 



j:oiq,iuj) = -uTJ2 J -^D{k,iiyn) (21) 



(Pk D{k + q, ivn + iuJn) 



(22) 



where: 



/cPp 
J^:^D{p, ivn) D{p + q, ivn 



(23) 



In this expression, D denotes the full propagator to order 1/A^'' (i.e including the self-energy Eq). Hence, Eq. H21() 
should be viewed as a self-consistent equation for Eq. We first consider this equation aX uj — q — defining a "gap" 




FIG. 2: Self-energy graphs in the expansion. 
A(r) — r — Y.{q = 0,iuj = 0) related to the correlation length f by A = A vanishes (^ diverges) at the T = Q 
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QCP. At the saddle point (N = oo) level A(T) = Eo(T = 0) - T,q{T). We show in Appendix El that, at this order, 
the self-consistent equation for the gap reads (see also Ref. 1251) : 

= i In + In + A,. In A,. - A,. (24) 

u 2 A + A,2 r(A + A,2) 9 9 * ' 

where Aq2 is a cutoff for g^, Aq2 = Aq2/(27rT) and A = A/(27rT). In the limit A < 1 and A < Aq2 F.q.^Q reduces 
to: 

^^A = -ilnA-AlnA (25) 

We dropped the index in A„2 as it is the only scale we use. In the limit u In A 1 the approximate solution of Eq. H25f) 
is^: 

ln(lnA/27rT) , , 

^--^ lnA/2.r 

Hence, A vanishes faster than T as T — > 0. The contributions of the marginally irrelevant coupling are smaller than 
temperature (or energy) itself. As we shall now sec, this also applies to the T — contributions to the Landau 
damping at the critical point, which is the main quantity of interest in this paper. 

We perform a calculation of the damping rate at the T — QCP up to order At the saddle-point level (l/A^^), 

the propagator at the QCP reads: D{q,iuj) = (|a;| +q'^)^^. The polarization bubble Uo{q,iLu) defined by can 
then be calculated exactly. Details are given in Appendix [CI Here we present the result: 

o 2tt . ■ N H + f, |w| + A ^, l^^l + g^X \uj\ + A 



g2 \uj\ + q^ + 2A-s \2\uj\ + q-^ J "\ 2\uj\ + q^ 



where s ~ \J (\io\ -V q^)"^ + 4g2A and Li2 (x) = — /q ln( 1 — y)/y . The above expression has the following asymptotics 
at uj/A 0, q^/A 0: 



g2 \uj\ 



4-21nV, 

A g2 



|a;| q^ u'^ 



2+ 21nV, ^^0, 

4 A \uj\ q^A 

TT^ \uj\ g2 g^A , , 

-2+ — -21nV. r^^O. (28) 



4 A 



\LiJ\ UJ'- 



We use the explicit expression (I27|) into H22|l and perform the frequency and momentum integrals numerically. The 
resulting self-energy Ei(g, jcj) is plotted in FigOl It is seen that the \/N corrections to the damping rate are less 
singular than lo at low-frequency. Alternatively, the corrections can be put in the form of an effective, scale-dependent 
exponent a(a;), defined by: a{ijj) = a;(9(lnx~^(iu;))/9u; with x~^(^i^) — \^\ — Si(g,«w = 0) -I- Si(g = ^ 0) the 
g = dynamical susceptibility at the QCP. This effective exponent is plotted on Fig. ^ it is seen to be only weakly 
dependent on frequency in the range where it is displayed. Hence, the logarithmic corrections in the marginal case 
can be mimicked by a power-law, but correspond to an effective exponent a > 1, in contrast to the experimental 
observation a < 1 for the two compounds mentioned above. 

Incidentally, we would like to comment on the data analysis made in Ref. 6 in order to support u/T scaling of the 
inelastic neutron scattering data in CeCus.gAuo.i. To this aim, we perform a similar analysis on our analytical result, 
by plotting T"x([j,T) versus uj/T (Fig.|3J). In this plot, the frequency-dependent susceptibility at finite temperature 
and g = is approximated by: x(cL',r)^^ ~ —iuj + A(T) — Ei(g = Q,uj) + Ei(g = 0,a; = 0) (this approximation 
does not include the temperature dependence of the Bose functions in the calculations which would be required to 
obtain the correct damping coefficient, as is required to obtain (|13|) ). Such a scaling plot is attempted for a varying 
range of a, and the optimal value of the exponent is chosen such as to yield the best coUapseS, as measured by the 
standard deviation plotted in Fig. It is seen that an excellent collapse is obtained with a ~ 1.09. We emphasize 
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0.3 




FIG. 3: Self-energy at the d = 2, z = 2 QCP, up to order as a function of (imaginary) frequency to. The calculation is for 

u/{2-Kf = 5, and LO is measured in units of the cutoff A. 



1.1 




1 I ' ' ' ' ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 



FIG. 4: Effective exponent a{LLi), as defined in the text, as a function of (imaginary) frequency uj. The parameters are as in 
FigEl 

that Lo/T scaling, in a strict asymptotic sense, does apply here, but with a trivial Gaussian exponent and scaling 
function: Tx{o->,T) = i/uj {uj — lo/T). The logarithmic corrections characteristic of the marginal case apparently 
mimic non-trivial scaling properties over quite an extended range of uj/T. This should serve as a warning for the 
interpretation of the experimental results. 

We end this section by considering the local approximation to the damped (fA theory in the marginal case, up to 
order In a previous paper, Motome and two of us^° have investigated the EDMFT approximation to the critical 

behavior of various models, and shown that this local approximation is quite satisfactory above the upper critical 
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FIG. 5: Scaling plot for r"Imx(w, T, g = 0) versus u/T. The temperatures displayed are: T = .03 — .9A. The choice a = 1.09 
provides a good collapse of the data over this frequency range. Different symbols correspond to different values of T. 




FIG. 6; Standard deviation criog(Tax)i rneasuring the quality of the lo/T collapse-, versus a. 



dimension (i.e here for d + z > 4), when the quartic coupling is irrelevant. In this approach, all skeleton graphs are 
taken to be local. For arbitrary z and d = 2, the local propagator at the T — QCP reads: 

f qdq 1 1 , /a+ \uj\i\ , , 

Performing the Fourier transform, this yields the asymptotic behavior at long (imaginary) time: 

D{t) = / 7—:;- In V e"'"^ « , / dxlnx"- cosx = (30) 

Loo 27r47r \^ \^\\ J {"2^^^ Jo ^7rz\T\ ^ ^ 

The 1/N correction to the spin-fluctuation self-energy is: 

St(r) = ^i^(r)3«-^^-^^ (31) 
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which corresponds to the following low-frequency behavior (with Ai a cutoff of order 1): 

1 3 

TV (4^(" 2 + 7 - li^^i + 1^ 1^1 V' (32) 

For real frequencies H31II corresponds to corrections to the Landau damping of the form ImSi(aj + jO+) oc uj'^Sign{uj). 
This is much smaller (by a factor of order oj, up to logarithms) than the corrections obtained from the direct calculation 
in d = 2 detailed above (Fig. (SJ- Hence, we conclude that the local (EDMFT) approximation is not very reliable at 
the upper critical dimension (marginal case). As we shall see in the following section, a numerical solution of the 
EDMFT equations for the Kondo lattice using an approximate "impurity solver" based on the large-N expansion 
leads to a similar conclusion: only subdominant corrections to the Landau damping rate are generated instead of the 
anomalous (power-law like) enhancement observed experimentally. 



IV. "LOCALLY" CRITICAL POINT: SELF-CONSISTENT LARGE-N SOLUTION OF THE EDMFT 

EQUATIONS 

A. EDMFT of the Kondo lattice and mapping onto a local spin-fluctuation model 

It has been recently arguedi^ii^ that the understanding of non-Fermi liquid behavior in heavy-fermion compounds 
close to a QCP requires a formalism in which electronic degrees of freedom and spin fluctuations can be treated on 
the same footing, at least as a starting point. This is a non-trivial task, which is however made easier in the context 
of the extended dynamical mean-field theory (EDMFT) ^^i^^i^''' of the Kondo lattice model, which we consider here. 

Let us consider the Kondo lattice Hamiltonian with an explicit exchange coupling between localized spins, chosen 
to be Ising-like: 

H = ^ ekclcr^ka + Si.c\^T„„,Cia' + ^ hj S^Sj (33) 

ka iaa' i<j 

In this expression, Jk is the (antiferromagnetic) Kondo coupling, and r stands for the Pauli matrices. EDMFT maps 
this model onto a local impurity model with effective action: 

S ^ ^ [dr [dT'Y,[cUT)g,\r-r')c„ir') + S^ir)xo'ir^r')S^T')] 
Jo Jo 

+ /drJK VS^.4f,.,c,, (34) 

Both Qq and xo are effective "Weiss fields" which must be determined in such a way that the following self-consistency 
conditions hold: 

/P (^) 
•^^^ FTT^ ^ (35) 

Xiiv,,) = <S\iv^)S'\~iv,,)>s^ f de /"'^f (36) 

J M[lVn)-€ 

In this expression, Sc(ia;„) and M(iun) are respectively the fermionic and local spin self-energies calculated within 
the impurity problem (|34ll . i.e.: 

G~^{iLUn) = Ga^iiuJn) - T,c{iLL!n) (37) 
X'^il^n) = -XQ\i'^n)+M{Wn) (38) 

The densities of states po and pi associated respectively with the conduction band and to the spin interactions read: 
Po(e) = J2k^(^ ~ ^k), Pi{^) — '^Zq^i^ ~ ^q)- The self-consistency conditions ()35I36|I simply express that the impu- 
rity model Green's functions and self-energies should coincide with their lattice counterpart (assuming momentum- 
independence of the self-energies). They can be recast into a more compact form, relating directly the local correlation 
functions to the two Weiss fields appearing in the effective action (|^ : 

G{iLo,,) = [ de ^ ^ (39) 

Xi^-n) = [ de P'^'l (40) 

J X (I'^n) +Xq {^'^nj - e 
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The key question at this stage is whether the electronic degrees of freedom can be integrated out, particularly near 
the QCP. The local EDMFT framework allows to handle this issue in a more controlled manner. Indeed, electronic 
degrees of freedom can be integrated outf from the impurity model (|34|l using an expansion in the spin-flip part 
of the Kondo coupling a la Anderson- Yuval-HammannS&Si. This can be conveniently performed, for example, using 
bosonization methodsi2i22i22. This leads to the following action for the spin degrees of freedom: 

S= dr dr'S'iT) [-XoHr - t') + K{t - r')] S'{t') + / dr T^" (41) 
Jo Jo Jo 

In this expression, K(t — r') ~ l/(''" — t'Y' is a long-range interaction induced by the electronic modes. It corresponds 
to the bare Landau damping {K{ivn) — K.\vn\ at low- frequency) . The coupling in front of the spin-flip term is 
proportional to the Kondo coupling: T oc tqJk (with tq a short-time cutoff). Hence, the problem has been mapped 
onto an quantum Ising model in a transverse field with long range interactions^^^iSi^. 

B. Numerical solution based on the large-N approximation 

This mapping has been used in order to perform analytical^^ and numerical-^- studies of the behavior at the QCP 
in mean-field models with random exchange couplings lij. It was demonstrated that this case is formally similar— 
to the Hertz-Moriya-Millis theory of an antiferromagnetic QCP in d — i. Recently, Grempel and Si have analyzed 
the EDMFT theory of the QCP in the two-dimensional case using this mappingiS. Their claim is that anomalous 
Landau damping (oc oj" , a < 1) and w/T scaling can be demonstrated in this case. Here, we present a solution of this 
problem in which the local impurity problem is solved using a large-N method. In contrast to the claims of Ref. il^ 
we find that only subdominant corrections to the Landau damping are generated. 

Following Ref. HtL we deal with the local action ()41|l for the transverse field Ising model by extending it to an 
iV-component rotor model n: 



1 



[drfi) 



+ A(T)(n2 -iV) 



i /dr /dr' n{T) [-Xo\r - r') + K{t - r')] n(r') (42) 



"'0 



In this expression, g plays a role to similar to that of the transverse field in the original Ising model (disordering 
term). We define a self-energy correction for the local spin susceptibility xi^) = {S'^{0)S^{t)): 

X{ivn) = ^-ri (43) 

Let us first consider the problem at the saddle-point level [N — oo), for a given Weiss function xo- The self-energy is 
then frequency- independent (as noted by Grempel and Sii^) and is given by the saddle-point value of the Lagrange 
multiplier: 

S]o(«i^n) = -A (44) 

such that: 

Xoo{t = 0) = 1 (45) 
In this expression, Xoo is the local susceptibility at the saddle-point {N — oo) level: 

^^^^''"^ " 2/ ^T^c / -!(■ (4^) 

Kl9 + K{lVn) - Xo yi-^n) + A 

Both A and Xoo depend on the Weiss function xo- We approximate the self-energy by its first two terms in the 1/N 
expansion: 

S(^^/„) ~ So(iz/„) + Si(i;/„)-Si(iO) (47) 
Si(r) = -^r(T)xoo(T) (48) 

with T{iur,) = — ^ and H(r) = [xoo(t)]' (49) 
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The EDMFT equations are thus solved by iterating (numerically) the following procedure. For a given Xoi'i'^n), the 
saddle-point quantities A = — Eq and Xoo are computed by solving Eqs. (|45I46|I . Then, the self-energy is computed 
from H47I48I49|I . This is inserted into the EDMFT self-consistency condition in order to obtain the local susceptibility: 



Xim.) = I Ae ^j,,.'''^ , (50) 

From this, an updated value of the Weiss function is obtained as: Xo^i'i'^n) — + K{ivn) — S(w„) — X~^(*i^n)j 
and the procedure is iterated until convergence is reached. 

The result of this numerical calculation for the (local) spin self-energy is displayed in Fig.[71 The behavior of the self- 
energy at low imaginary frequency is well fitted by S]i(iaj) ^ In |a;| (inset) when one tunes the parameters to sit at 
criticality (Fig. (SJ. This corresponds to subdominant corrections to the Landau damping ImE(a; + oc a;^sign(a;). 
These findings are in agreement with the expansion of the EDMFT equations for the 0'* model at order 1 /N discussed 
at the end of the previous section, and provide an independent check of this result. The method followed in the 
present section has been to use the large-N method as an approximate "impurity solver" for the local problem. It 
can be easily checked that if, instead, the equations used here are all expanded up to order 1/iV, the equations of the 
local approximation presented in the previous section are recovered. 




FIG. 7: Self-energy E(zaj) close to the critical point (at inverse temperature (3 = 500 and g = gc = 0.223), as a function of 
imaginary frequency. For this calculation, we have set N = 2, and used a flat density of states pi{e) of half- width equal to 
unity. Inset: plot of ti;~^E(iaj), showing the logarithmic behavior of this quantity at small frequencies. 



V. CONCLUSION 



We analyzed two of the approaches to a critical theory of the antiferromagnetic metal to paramagnetic metal in 
the CeCu6_2:Au^ and YbRh2Si2 materials from the perspective of the antiferromagnetic spin fluctuations. 

In the context of the model of Rosch et. ali^ we analyzed the corrections to mean field theory due to the marginally 
irrelevant fourth order coupling. We find that while there are logarithmic corrections to mean field theory which 
vanish as l/Log^|cij| the coefficient of this corrections is such that it can only be mimicked by a decreased Landau 
damping, in clear disagreement with experimental observations. In our view, this rules out this model, in spite of the 
fact that it can account nicely for many of the electronic properties of this systems. 

We also analyze the "local " version of this model using the large N expansion. Again in this case, we find that it is 
not possible to account for the anomalous frequency dependence of the damping coefficient. This is in disagreement 
with referenceiSi but in agreement with a recent EDMFT study of the Anderson lattice modeli^ 

The anomalous damping has been recently observed in NMR measurements in both^SiS YbRh2Si2 and 
CeCug-xAux. While it has different temperature dependencies in each system, it is clearly enhanced above the 
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0.25 



^ 0.15 - 




0.05 



FIG. 8: Plot of the mass E(iO) as a function of temperature, demonstrating the existence of a continuous transition at T = 
within our approximate solution of the EDMFT equations. From bottom to top: g = 0.1 (ordered state), g — 0.223 (critical at 
r = 0), (? = 0.3 (disordered case). 

fermi liquid temperature dependence. Possible resolution of this puzzle may require a better treatment of the effects 
of disorder or better treatments of models which retain the fermionic nature of the spin fluctuations as proposed by 
Coleman et. al,^^ and Si et. alH 
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APPENDIX A: COMPUTATION OF RG FLOW 

The computations leading to the RG flow © are standard. Here we only provide some further details of the 
computation of the number C2, and of the constant c required to obtain the universal logarithmic correction in 
These computations involve the |a;| term in the propagator in an essential and novel manner. 

Repeating the standard field-theoretic derivation for the two- loop graph in Fig ^ we have: 



hiq^ii^ = 0) 



<{N + 2) 
18 



duji duj2 



d'^ 



92 



((q + qi + q2)^ + l^^i + ^^2!) 



(Al) 



I2 is to be evaluated for d < 1, where it is convergent, and the result analytically continued to near c? = 2, where it is 
expected have the small e expansion shown on the r.h.s. The coefficient of the pole fixes the value of the constant C2. 

Note that we are evaluating I2 in an external momentum q, and in d = 2 the pole on the r.h.s. corresponds to 
the appearance of q^lnq term. In contrast, evaluation of 12(0,10;) (presented below) shows that no pole appears, and 
hence there is no \uj\ In |a;| term in d = 2. This absence of such a pole is, of course, the reason for the absence of an 
independent renormalization constant for the damping term noted below lO, and for the u dependence of I2 noted 
above It is also responsible for the exponent identity 2 = 2 — ?7for(i<2. 
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We now turn to the evaluation of I2 in (jAl|) in an external momentum q. First, we evaluate the integral over 52 by 
the standard Feynman parameter method 



ul{N + 2)T{2- d/2) 
18(47r)'i/2 





duJ2 


f d% f 


J -00^ 


'2^ J 


C^ttY Jo 




1 





dx- 



1 



[(q + qi)2a;(l - x) + \uj2\{l - x) + \uj^ + LJ2\xf"''^ ' 
Similarly, performing the integral over qi with a Feynman parameter y we obtain 

ug(7V + 2)r(3-d) f°° duJidLj2 



(A2) 



/2('7,0) 



18(47r)° 



2tt 2tt 



dx[x{l - I 

1 



[q^x{l ~ x)y{l -y) + \lui + lu2\x{1 - x){l -~ y) + \uji\y{l - x) + \LU2\xy]^+' 
Now we perform the integral over ui and UJ2 by using the useful formula 

1 _ 4A^-^{B + C + D)T{a -2) 



duJidiLj2 



[A + B\uji\ + C\uj2\ +D\oji+ 1^21]" {B + C){C + D){B + D)r(<T) ' 



(A3) 



(A4) 



The formula l|A4p is derived by explicitly performing the integrals over lji_2 over different regions in the 0^1,2 plane 
delineated by changes in signs of uji, uj2, and uji + uj2- Note that (|A4p is to be evaluated at cr = 1 + e. Consequently, 
we obtain from (|A4|) the factor r{a ~2) = r( — 1 + e) = — l/e + 0{e'^), which gives us the requisite pole in e appearing 
on the r.h.s. of (|A1() . So we may safely set d = 2 in the remaining terms in (|A3() . In this manner, combining HA1(I . 
((X3|) and HA4|I we obtain 



C2 



(iV- 



36 



dx 



\ {l-y) [x{l - x){l -y) + y(l - x) + xy] {N + 2)(12 - -k^) 
ay- - 



[x{l - y) + y][{l - x){l ~ y) + y] 



144 



(A5) 



Let us now present a few details of the evaluation of /2(0, ioj). In this case, (|A3|) is replaced by 



/2(0,ZW) 



ul{N + 2)T{3-d) 
18(47r)'* 



/ dx[xil-x)Y/^ 
Zn Zn ./n 



1 



dyy 



e/2 



[\uj + UJI+ uj2\x{l ~ x){l - y) + \uji\y[l - x) + \L02\xyY+^ ' 
Now, instead of ljA4p . we need the following integral 

1 _A\uj\'^-''ABCT{a -2) 



(A6) 



du!idui2 



[A\uJi\ +B\l02\ + C|W + +UJ2\] 



(A2 -B2)(A2 - C2) (B2 _^2)(52 _ (72) (^2 _ ^2) ((72 _ ^2) 



(A7) 



As below (jA4|l . we need to evaluate (|A7|l at a — 1 + e and pick out a possible pole in e. Indeed, a possible pole does 
appear to be present in the r{a — 2) pre-factor in (|A7p . However, careful evaluation shows that that the residue of 
such a pole vanishes, and ljA7|l in fact has a smooth cr — > 1 limit: 



limjSZj = AABClul 

(T—*l 



In A 



InB 



+ 



InC 



(A2 -S2)(yl2 -C2) (52 -yl2)(S2 -C2) (C2 -^2)(C2 _52) 



(AS) 



The e — > limit of the remaining terms in ljA6p is straightforward, and this establishes the claimed absence of a 
|w| In |a;| term in /2(0, ii^) in d — 2. Evaluating the constant c in /2(0, iuj) — /2(0, 0) — cg'^\LL>\ we obtain: 



dx / dy 



x2(l-y) 



(l-2x) [2/2 -a;2(l -y)2 



■In 



(iV + 2)(67r2 ln2- 11C(3)) 
1728 



(A9) 



We performed integrations in Ea. ljA9|) by changing variables to x' — x, y' — x{l — y)/y and integrating first over x' 
and then over y' . 
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APPENDIX B: GAP EQUATION 

In this appendix we derive the gap equation Ea. (|^ . We need to evaluate: 

Eo{T) = -uTY I '^-^D{qM = ^^rVln , ^'^^ (Bl) 
£-^7 27r ' 47r ^ |cj„|+A + A,2 ^ ^ 

We consider a cutoff A^, which from now on we denote as A, and a frequency cutoff A^^ G 2TrnT which is to be 
taken to infinity. For this choice of the sharp cutoff A^^ one can show: 

SAT)= V m, + ^ =lnA±A + 21n^(^ + ^J^(^ + ^") (B2) 

-A.^:;<A. Kl + A + A,. A r(A)r(A + A + A„) 

where A = K/{2nT), = A^/{2TrT) and A ^ A/{2ttT). Taking the hmit A^ c» we have: 

(T) = In A±A + 2 h. - 2 A In A. (B3) 

^ ^ A r(A) ^ ' 

The finite T self energy is Eo(T) = ^wTS'i(T), and the zero temperature So(T = 0) is: 



So(r = o) = ^wr^i(r) 

Noticing that A(r) = Eo(r = 0) - Eo(7') we write Ea. ipijl . 



uA „ A , ,„ , 

.(In— -1) (B4) 



.0 (2^)^' A. 



APPENDIX C: POLARIZATION BUBBLE Do 

Here we show in detail how we compute Tlo{q, iuj) in d = 2 with z — 2. We use a finite cutoff Ag2 = A for in a 
momentum integration, and we use A^ — > oo as a frequency cutoff. We use the Feynman parameterization to find: 



Ilo{q,iu) 



(Pp dv 
(27r)2 2^ 



Z3(p, iv)D{p + q,iiy + ilj) 



4vr 7o 



dx 



+ OC 



dv 
2^ 



(l{x,q,v,iuj) ^ - [A + l{x,q,i',iuj)] ^) 



(CI) 



where l{x, q, v, iuj) — x{l — x)q^ + x\v\ + (1 — x)\v + u)\ and D ^ = + 9^ at the QCP. Integrating out frequency we 
get: 



U,{qM = ^,vj^ dx- 



Ax 



2x 



{In [{l-x){q^x+\u\ 



In 



— {q'-\Lu\ + s)-x][ q'x + -q' + s) 



(C2) 



where s — a/ (|ti;| + q^Y + Aq^A. It is easy to show: 



dx—— In ax + 6=2- 2 In - 21n6 - L12 — — + Lia — - 

1 — 2x a b a + 2o a + 2b 



(C3) 



for a + 6 > and 5 > 0. From Ea.l(n2)) and Ea. (|n3)l . after simple algebra, Eq.||22I) follows. 
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